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Abstract 

Tension leg platforms are used for oil exploitation in deep 
water. The accurate and reliable responses are needed for 
optimum design and control of the structure. Free vibration 
analysis gives us a deep view of structural response and the 
judgment on coupling effects between degrees of freedom. 
In this paper the nonlinear dynamic analysis of TLP is 
carried out. The power spectral densities (PSDs) of responses 
are calculated from nonlinear responses. A computer 
program has been developed capable of solving differential 
equation of motion of TLP under various types of excitations. 
The focus of the paper is on the conceptual discussion of the 
structural behavior and responses of the structure related to 
initial arbitrary conditions. As an example a case study is 
investigated and numerical results are discussed. 
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Introduction 

It is obvious that there is an increasing demand for oil 
exploitation in deep water. As the water depth 
increase, the environment will be severer and therefore 
some innovative structures are required for economic 
production of gas and petroleum in deep water. There 
are several systems that have been developed for deep 
water. Because of high stresses if fixed structured used 
in deep water, deep water systems will be compliant 
and therefore these structures must be designed 
dynamically. One of the important types of these 
structures is TLP (Fig 1) 

Many studies have been carried out to understand the 
structural behavior of TLP and determine the effect of 
several parameters on dynamic response and average 
life time of the structure (Faltinsen et al., 1982; 
Angelides, et al., 1982; Teigen 1983; Ahmad, 1996; Jain, 
1997; Chandrasekaran and Jain, 2002; Tabeshpour et al., 
2006a). The tether system is a critical and basic 
component of the TLP. The most important point in 
the design of TLP is the pretension of the legs. The 


pretension causes that the platform behaves like a stiff 
structure with respect to the vertical degrees of 
freedom (heave, pitch and roll), whereas with respect 
to the horizontal degrees of freedom (surge, sway and 
yaw) it behaves as a floating structure. Therefore the 
periods of the vertical degrees of freedom are lower 
than that of the others. 

Another important task is to investigate the effects of 
radiation and scattering on the hull and tendon 
responses. An analytical solution to surge motion of 
TLP has been proposed and demonstrated (Lee et al., 
1999), in which the surge motion of a platform with 
pre-tensioned tethers was calculated. In that study, 
however, the elasticity of tethers was only implied and 
the motion of tethers was also simplified as on-line 
rigid-body motion proportional to the top platform. 
Thus, both the material property and the mechanical 
behavior for the tether incorporated in the TLP system 
were ignored. 

An important point in that study is linearization of the 
surge motion. But it is obvious that the structural 
behavior in the surge motion is highly nonlinear 
because of large deformation of TLP in the surge 
motion degree of freedom (geometric nonlinearity) 
and nonlinear drag forces of Morison equation. 
Therefore the obtained solution is not suitable for the 
actual engineering application. For heave degree of 
freedom the structural behavior is linear, because there 
is no geometric nonlinearity in the heave motion 
degree of freedom and drag forces on legs have no 
vertical component. Similarly, an analytical heave 
vibration of TLP with radiation and scattering effects 
for damped systems has been presented (Tabeshpour 
et al., 2006b). A similar method is presented for 
hydrodynamic pitch response of the structure 
(Tabeshpour et al., 2006c). The modified Euler method 
(Karman and Biot, 1940) presented herein is a simple 
numerical procedure which can be effectively used for 
the analysis of the dynamic response of structures in 
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the time domain. It has been shown that the modified 
Euler method is conditionally stable (Hahn, 1990). 

The application of the modified Euler method made 
herein shows that it is efficient and easy to use, and 
that it can be employed to obtain accurate solutions to 
a wide variety of structural dynamics problems. 
Simplicity is one of the distinguishing features of the 
method. Because the modified Euler method is 
conditionally stable, it may be inefficient for the 
analysis by means of direct integration of the response 
of a multi degree-of-freedom system with a very short 
but the highest natural period of vibration. However, 
the method is explicit, and particularly suitable for the 
analysis of non-linear systems. The modified Euler 
method has been successfully used in the analysis for 
the dynamic response of wave-excited offshore 
structures (Sanghvi, 1990). 

The effect of added mass fluctuation on the heave 
response of tension leg platform has been investigated 
by using perturbation method both for discrete and 
continuous models (Tabeshpour et al., 2006d). 

Liu et al. described an analysis of the non-linear effects 
and identification of non-linear pitch motion on 
tension leg platforms. The purpose of their paper was 
to accurately identify pitch motion on the tension leg 
platform and to interpret the non-linear effects using 
statistical methods, the NARMAX methodology, and 
the higher order frequency response functions (Liu et 
al., 2004). 

A robust stochastic design framework has been 
discussed for design of mass dampers by Taflanidis et 
al. (2009) was focused on applications for the 
mitigation of the coupled heave and pitch response of 
Tension Leg Platforms under stochastic sea excitation. 

Tabeshpour et al. (2010) investigated design and effect 
of tuned mass damper on response of tension leg 
platform under wind and wave forces. Efficiency of 
numerical methods is a priority to dynamic nonlinear 
analyses. It is possible to use numerical methods such 
as beta-Newmark in order to investigate the structural 
response behavior of the dynamic systems under 
random sea wave loads but because of necessity to 
analyze the offshore systems for extensive time to 
fatigue study, it is important to use simple stable 
methods for numerical integration. The key point of 
suitability of MEM to solve the TLP system is that the 
maximum frequency of the system is about 0.5 Hz. 
Tabeshpour et al. (2006e) has investigated the 
efficiency of the modified Euler method (MEM) for 


nonlinear dynamic analysis of tension leg platforms as 
a response dependent system. 

Convergence and stability analyses of MEM were 
presented by Tabeshpour et al. (2012) Modified Euler 
method (MEM) has been used for solving the 
differential equation of. Surge motion analysis of TLP 
under linear wave via perturbation method has been 
presented by Tabeshpour and Shoghi (2011) who has 
also been presented limitation between linear and 
nonlinear models for surge motion of TLP. 

In this paper an integrated response analysis of TLP 
under arbitrary initial values is presented giving a 
deep view of results. Displacements and strains of 
tendons are investigated and discussed in this paper 
under initial arbitrary independent values. Modified 
Euler method is used to solve nonlinear differential 
equation of motion. Phase plane used in here is a 
helpful tool to discuss the stability and convergence of 
the response. 

Structural Modelling of TLP 

Mass Matrix 

Structural mass is assumed to be lumped at each 
degree of freedom. Hence, it is constant and diagonal 
in nature. The added mass, Ms, due to the water 
surrounding the structural members and arising from 
the modified Morrison equation, is considered up to 
the mean sea level (MSL) only. The fluctuating 
component of added mass due to the variable 
submergence of the structure in water is considered in 
the force vector depending upon whether the sea 
surface elevation is above (or) below the MSL. The 
mass matrix of TLP is 
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mass moment of inertia about the x axis = Mr, 2 , M pp is 
the total mass moment of inertia of the y axis = M r 2 , 
M YY is the total mass moment of inertia of the z axis = 
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Mr. 2 , Tx is the radius of gyration of the x axis, r y is the 
radius of gyration of the y axis, and Tz is the radius of 
gyration of the z axis. The added mass terms are: 

M aSS =M aWW =M aHH = 0.25nD 2 (C m - l)p dl (2) 

M ass = j* dM a ss (3) 

lenght 

M MS , M clRW and M aRn are the added mass moment of 
inertia in the roll degree of freedom due to 
hydrodynamic force in the surge, sway and heave 
direction, respectively. M aPS , M rit , w , and M aP „ are the 
added mass moment of inertia in the pitch degree of 
freedom due to hydrodynamic force in the surge, sway 
and heave direction, respectively. The presence of off 
diagonal terms in the mass matrix indicates a 
contribution to the added mass due to the 
hydrodynamic loading. The loading will be attracted 
only in the surge, heave and pitch degrees of freedom 
due to the unidirectional wave acting in the surge 
direction on a symmetric configuration of the platform 
about the x and z axes). 



FIG. 1 TLP CONFIGURATION 
Stiffness Matrix of the TLP 

The coefficients, Kab, of the stiffness matrix of the 
triangular TLP are derived as the reaction in the 
degree of freedom A due to unit displacement in the 
degree of freedom B, keeping all other degrees of 
freedom restrained. The coefficients of the stiffness 
matrix have nonlinear terms due to the cosine, sine, 
square root and squared terms of the displacements. 
Furthermore, changes in tendon tension due to the 


motion of the TLP in different degrees of freedom 
make the stiffness matrix response-dependent. The 
stiffness matrix [K] of a TLP is: 
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S, W, H, R, P and Y denote surge, sway, heave, roll, 
pitch and yaw respectively. These elements have been 
investigated by many researchers [2-7] . 

In the stiffness matrix, the presence of off-diagonal 
terms reflects the coupling effect between the various 
degrees of freedom and the coefficients depend on the 
change in the tension of the tendons, which affects the 
buoyancy of the system. Hence, the [X] is not constant 
for all time instants but the coefficients are replaced by 
a new value computed at each time instant depending 
upon the response value at that time instant. 

Damping Matrix, [C] 

Assuming [C] to be proportional to [X] and [M], the 
elements of [C] are determined by the equation given 
below: 

C = aM + (3X (5) 

a and (3 are constant. This matrix is calculated based 
on the initial values of [X] and [M] only. There are two 
approaches to define a and (3 . The first method to 
calculate them in a manner that a predefined damping 
ratio E, = C / 2 M to is achieved in which to is angular 
frequency of the system. The second method is to 
consider an appropriative value for them such as 2%- 
5%. 

Equation of Motion 

The free vibration of equation of motion of the TLP is 
given as: 

[M ] [X } + [C ] {X } + [K ] [X } = {F (0} (6-a) 

{F(D} = {|0.5 Plv C d Dx |i| + (0.257tD 2 [C,„ -l]p„,i')} (6-b) 

where [M], [C] and [X] are the matrices of mass, 
damping and stiffness respectively, {X } , {X } and 
[X } are the structural displacement, velocity and 
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acceleration vector respectively and (F(t)} is the 
excitation force vector. 

Case Study 

A TLP in 800m deep water has been chosen for the 
numerical analysis. The characteristics of the TLP are 
shown in table 1. Free vibration analysis is carried out 
using a MATLAB based code developed for nonlinear 
dynamic analysis of TLP. 


TABLE 1 TLP PROPERTIES 


Platform height 

80 m 

Platform weight 

1.3x10 s IV 

Total pretension 

0.8x10 s IV 

Draft 

45.28 m 

Comer column diameter 

18 m 

Center to center column spacing 

90 m 

Pontoon length 

90-18 m 

Pontoon diameter 

12 m 

Center of gravity above the 
base line 

10 m 

Tether length 

800-45.28 m 

Diameter of tether pipe 

1.2 m 

Young s modules of tether 
material 

2.019x10“ JV / m 2 

Number of columns 

4 

Service life 

20 years 

Radii of gyration (R x ,R y ,R v ) 

55 m, 55 m and 65 
m 

Drag coefficient, C D 

1.0 

Inertia coefficient, C 

1.8 

Mass density of water 

1025 kg / m 3 


Draft is calculated as follows: 


(63.39s) is lower than that of sway (76.67s). Table 3 
shows TLP periods and frequencies for all degrees of 
freedom. 


TABLE 2 TLP PROPERTIES 


Motion 

Initial Value 

surge 

30 m 

sway 

20 m 

Heave 

0.5 m 

Roll 

0.01 rad 

Pitch 

0.01 rad 

Yaw 

0.01 rad 


TABLE 3 TLP PERIODS AND FREQUENCIES 


Motion 

Period (sec) 

Frequency (rad/s) 

surge 

63.39 

0.1 

sway 

76.67 

0.08 

Heave 

1.98 

3.17 

Roll 

1.74 

3.61 

Pitch 

1.74 

3.61 

Yaw 

80.57 

0.078 


Fig 2 shows the nonlinear relation between the forces 
in surge direction and surge displacement. The system 
is hardening and geometry nonlinearity is important 
to increase the stiffness in large displacements. Surge 
stiffness is as follow: 
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where W is the platform weight, T is pretension of 
tendons, p is water density, D c is diameter of column, 
D is pontoon diameter, S is distance between 
pontoons. 

Initial values assigned for free vibration of TLP are 
shown in table 2. Because the stiffness is related to 
displacement and initial displacement of surge motion 
(30m) is more than sway (20m), therefore surge 
stiffness is more than sway stiffness and surge period 


It means that the stiffness is a parabolic function. 



Surge (m) 

FIG. 2 RELATION BETWEEN SURGE FORCE AND SURGE 
DISPLACEMENT 
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Surge stiffness has two parts: 


Weakly nonlinear part: 


nT n 


( 2 A 

l + x 


21 - 


This part is approximately linear. 
• Strongly nonlinear part: 


nTrAE x 2 
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, x 

1 H y 

2 r- 
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This part is parabolic. Figs 3 and 4 show weekly and 
strongly nonlinear parts of surge stiffness respectively. 

Fig 5 shows the surge stiffness versus surge motion. It 
is seen that in approximately 3% relative displacement 
(surge=27m) stiffness is becoming more than two times 
comparing zero displacement. 



(a) 



(b) 

FIG. 3 WEEKLY NONLINEAR PART OF SURGE STIFFNESS 



Surge (m) 

FIG. 4 STRONGLY NONLINEAR PART OF SURGE STIFFNESS 



FIG. 5 SURGE STIFFNESS VERSUS SURGE MOTION 

Fig 6 shows the variation of surge stiffness versus time. 
Comparing weak and strong parts of stiffness shows 
the importance of strong part in large displacements. 
The total varying stiffness in time will reduce to a 
minimum value for each cycle with two rounds. 



■■■■ weakly nonlinear response dependent 

strongly nonlinear response dependent 

— total 



2 



Time 


FIG. 6 VARIATION OF SURGE STIFFNESS VERSUS TIME 

Fig 7 shows time histories of responses. It is seen that 
all degrees of freedom start from initial dictated value 
to the zero point, but there is a short time of instability 
in heave motion because of the absence of unison and 


49 


www.daoe-journal.org 


Development and Applications of Oceanic Engineering (DAOE) Volume 2 Issue 2, May 2013 


compatibility among the arbitrary initial values. 
However the response converges to the zero 
equilibrium point. 

Heave response of left hand comer of TLP is shown in 
Fig 8. Transient response at the first three seconds 
shows a temporary instability of the response due to 
the exclusion of compatibility among initial values. 
However the response converges because the system is 
hardening. 




Time (sec) Time (sec) 





There is an initial strain in tendons because of 
pretention stress. Therefore when tendons are under 
tension in nonzero condition, the strain can be divided 
into two parts: initial strain and fluctuating strain. 
Fluctuating part of strain in right hand side tendon is 
shown in Fig 9. Similarly Fig 10 shows total strain in 
left hand side tendon. 



Time (sec) 

FIG. 9 FLUCTUATING PART OF STRAIN IN RIGHT HAND SIDE 
TENDON 



FIG. 10 TOTAL STRAIN IN LEFT HAND SIDE TENDON 


FIG. 7 TIME HISTORIES OF RESPONSES 



FIG. 8 HEAVE RESPONSE OF LEFT HAND CORNER 


Trajectory of TLP motion in space gives a global view 
of the dynamic response and the path of the system 
wending from the start point to zero point (Fig 11). 

A horizontal view of TLP motion is shown in Fig 12. 
the convergence from the start point (30, 20) to the zero 
equilibrium point can be observed clearly. 

Phase planes of all degrees of freedom show 
converging to unique equilibrium stable point. There is 
a special case in heave motion. Firstly the initial point 
is approaching from the initial point (1) to point 5 that 
has the largest distance from the initial point. This 
event repeats for some other cycles which happen 
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because the initial values are arbitrary and 
independent to each other. Then the point is 
approaching to the zero equilibrium point on the spiral 
path. 




FIG. 11 TRAJECTORY OF TLP MOTION IN SPACE 



Surge (m) 

FIG. 12 HORIZONTAL VIEW OF TLP MOTION 


Power spectral densities of sway and pith responses 
are shown in Figs 14 and 15. Because of coupling, there 
are two peaks in some degrees of freedom. Peaks are 
occurred at the related frequency of each degree of 
freedom. 



Surge (m) Ftoll (m ) 



FIG. 13 PHASE PLANES 



Angular Frequency (rad/sec) 

FIG. 14 POWER SPECTRAL DENSITY OF SWAY 



Angular Frequency (rad/sec) 

FIG. 15 POWER SPECTRAL DENSITY OF PITCH 
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Conclusions 

A conceptual discussion on the results of nonlinear 
dynamic analysis of TLP under arbitrary initial values 
has been presented. A special case has been examined 
at the initial transition part of heave motion diverging 
from equilibrium point. Because the system is 
hardening, this incompatible initial value is free from 
diverging problem. Modified Euler method suitable 
method for nonlinear dynamic analysis of response 
dependent systems was used to solve nonlinear 
differential equation of motion. The power spectral 
densities (PSDs) of responses were calculated from 
nonlinear responses. 
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